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ABSTRACT 

High velocity outflows from quasars are revealed by the absorption signatures they 
produce in the spectrum of the quasar. Clues on the nature and origin of these flows 
are important for our understanding of the dynamics of gas in the central regions of the 
Active Galactic Nucleus (AGNs) but also of the metal enrichment of the intergalactic 
space. Line radiation pressure has often been suggested to be an important process 
in driving these outflows, however no convincing evidence has been given so far. Here 
we report observation of a highly structured flow, toward Q 1511+091, where the 
velocity separations between distinct components are similar to O vi, N v and C iv 
doublet splittings with some of the profiles matching perfectly. This strongly favors 
the idea that the absorbing clumps originate at similar physical location and are 
driven by radiative acceleration due to resonance lines. The complex absorption can 
be understood if the flow is highly collimated so that the different optically thick clouds 
are aligned and cover the same region of the background source. One component shows 
saturated H i Lyman series lines together with absorptions from excited levels from 
C II and Si 11 but covers only 40% of the source of continuum. The fact that clouds 
cover only part of the small continuum source implies that the flow is located very 
close to it. 

Key words: Quasars: absorption lines - Quasars: individual: Q 1511+091 



1 INTRODUCTION 

Comprehensive understanding of the nature and origin of 
gas flows in the very central regions of AGNs is elusive de- 
spite decades of intensive investigation. Presence of outflows 
with very large velocities (~20000 km s^^) is revealed by the 
broad as well as narrow absorption lines (BALs and NALs) 
observed in the spectrum of AGNs. It has been suggested for 
long that the momentum transfered through resonance ab- 
sorption lines could play an important role in driving these 
flows as it does in the case of O star winds (Castor, Abbott 
& Klein 1975; Arav, Li, & Begelman 1994). However, signa- 
tures of acceleration, though claimed to be noted in nearby 
Seyfert galaxy NGC 3515 (Hutchings et al. 2001), are not 
generic. 

* Based on observations carried out at the European Southern 
Observatory (ESO) under prog. ID No. 65.P-0038 and 65.O-0063 
with the UVES spectrograph on the Very Large Telescope (VLT) 
at Cerro Paranal Observatory in Chile. 



A possible way to reveal line-driven acceleration is the pres- 
ence of absorption features with velocity separation similar 
to doublet and/or multiplet splittings (a consequence of the 
so-called "line-locking" effect). Such a structure in absorp- 
tion is achieved when there is a reduction in the flux that 
drives the far away gas element due to line absorption pro- 
duced by the gas closer to the QSO (Scargle 1973; Braun & 
Milgrom 1989). For higher efficiency, the gas from different 
elements of the flow must cover the same region of the back- 
ground source. Since the gas producing intrinsic absorption 
is known to cover the background source only partially (Pe- 
titjean et al. 1994; Hamann 1997; Barlow & Sargent 1997), 
line-locking can be achieve in this case if the flow is somehow 
collimated. This may be one of the reasons why tentative 
evidences for line-locking (Foltz et al. 1987; Srianand 2000; 
Srianand & Petitjean 2001; Vilkoviskiji & Irwin 2001) as well 
as double trough (the so-called Lyman-alpha ghost) seen in 
the mean profile of BALs (Weymann et al. 1991; Korista 
et al. 1993; Arav, Li & Begelman 1994) are not completely 
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Figure 1. The left-hand side panels show the fit to the eiTective continuum (continuum plus broad emission line) of QSO spectrum in the 
C IV (top) and Si IV (bottom) emission line regions. Right-hand side panels show the corresponding normalised fiuxes on a velocity scale 
with respect to z = 2.878. The y-axis in the left and right panels give the observed spectrum (without flux calibration) and normalised 
flux respectively. The continuous and dotted profiles are for the first (C IVA1548 and Si lvA1393) and second (C IVA1550 and Si lvA1402) 
members of the doublets respectively. The horizontal dashed lines in each panels give the residual flux expected for different covering 
factors (as indicated above each line). If the absorbing gas covers the continuum source only then the expected total covering factor of 
the background source is ~ 0.5 and ~ 0.70 for the C IV and Si IV lines respectively. These values are obtained using the relative flux 
ratios of emission lines and continuum at the position of the absorption line from low dispersion spectra of Sargent et al (1988). 



convincing. In this paper we present clear evidence for a flow 
that is highly structured through radiative pressure modu- 
lated by absorption in different resonance lines. 



2 OBSERVATIONS 

We have used the Ultra-violet and Visible Echelle Spectro- 
graph UVES (D'Odorico et al. 2000) mounted on the ESO 
KUEYEN 8.2 m telescope at the Paranal observatory on 
April 4-7 2000 to obtain a high-spectral resolution spectrum 
of the Zem = 2.878 and mv = 17 quasar Q 1511-1-091. Two 
standard dichroic settings have been used to observe with 
both arms of the spectrograph at the same time. The fi- 
nal wavelength range is 3260—10,000 A. The slit width was 
1 arcsec resulting in a resolution of ~45000. The total ex- 
posure time was 10 hours in seeing conditions better than 
0.8 arcsec full width at half maximum. The data were re- 
duced using the UVES pipeline, a set of procedures imple- 
mented in a dedicated context of MIDAS, the ESO data 
reduction package. The main characteristics of the pipeline 
is to perform a precise inter-order background subtraction 
for science frames and master flat-fields, and to allow for an 
optimal extraction of the object signal rejecting cosmic ray 
impacts and performing sky-subtraction at the same time. 
The reduction is checked step by step. Wavelengths were 



Table 1. Covering factor of the absorbing gas 



Component 




species 


fc 


a 


2.8186 


C IV 


0.10±0.04 


c 


2.8398 


C IV 


0.48±0.02 


d 


2.8523 


Ly 


0.30±0.01 






Ly series 


0.30-0.40 






C IV 


0.50±0.01 






Si IV 


0.30±0.05 






N v*" 


0.46±0.01 


e 


2.8607 


Ly 


^ 0.42 






Ly series 


0.40-0.45 






Al III 


0.43±0.03 






C IV 


^ 0.57 






Si IV 


^ 0.49 


f 


2.8641 


C IV 


0.77±0.01 






N V 


0.75±0.02 


J (blue) 


2.8908 


Ly a 


^ 0.15 






C IV 


^ 0.10 






Si IV 


0.28±0.12 






N V 


0.22±0.03 


j(red) 


2.8934 


Ly a 


^ 0.15 






C IV 


^ 0.10 






Si IV 


0.17±0.07 






N V 


0.18±0.02 



Line saturation is assumed 
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Figure 2. C IV absorption profile (top panel and solid line in other panels). The same profile, but shifted by different velocity separations, 
is overplotted as a dotted line in the other panels (respectively O VI, N V and C IV from top to bottom). Different velocity components 
that are separated by C IV (circle), N V (star) and O VI (filled square) doublet splittings are marked in the top most panel. The vertical 
dotted lines mark the position of various distinct component identified in the flow. 



corrected to vacuum-heliocentric values and individual ID 
spectra were combined together. This resulted in a S/N ratio 
per pixel of 25 at ~3600 A and 80 at ~6000 A. Typical er- 
rors in the wavelength calibration are less than ~0.5 km s~^. 



3 COVERING FACTOR 

The UVES spectrum of Q 151H-091 over the C iv and 
Si IV emission line wavelength ranges is shown in the left- 
hand side panels of Fig. ^ together with the fit to the effec- 
tive continuum (i.e., quasar continuum plus broad emission 
hues) we used to normalize the spectrum. Partial coverage 
of the absorbing gas implies some residual at the bottom 
of absorption lines even if they are heavily saturated. Par- 
tial coverage is apparent when the two resolved lines of a 
doublet, with different oscillator strengths, have the same 
apparent residuals without going to zero. This is obviously 
the case for most of the components in the C iv flow toward 
Q 1511+091. The two absorption lines of the C iv doublet 
corresponding to components a, c, d and f have exactly the 
same apparent depths over their profiles. It can be seen on 
right panels of Fig. ^ that when the spectrum is shifted by 
the velocity separation of the two lines of the C iv doublet, 
the profiles coincide nearly perfectly for most of the wave- 
length range. Using the method described in Srianand & 
Shankaranarayanan (1999) we estimated the covering factor 
for different species in individual components. The results 
are given in Table. |l| 

It is possible to say that different part of the flow cover the 



same background region. Indeed, it can be noted that when 
C IVA1550 of one component, say component 1, overlaps 
with C IVA1548 of another component, say component 2, 
the resultant optical depth, although smaller than 0.5 is not 
the sum of the optical depths of the two absorption lines. 
This is particularly striking for C lvA1548 of component 
e and C lvA1550 of component d or for C lvA1548 of 
component g and C lvA1550 of component e. This can be 
understood if one of the saturated components, say com- 
ponent 1, is shielding the other component, say component 
2 in such a way that the flux supposed to reach compo- 
nent 2 in C IVA1548 is absorbed completely by compo- 
nent 1 in C IVA1550. This explains why, although blending 
is apparent, the C lvA1548 profile matches so closely the 
same profile shifted by the velocity splitting between the two 
C IV lines (see Fig. 0). This shows that at least part of the 
flow is collimated and covers the same part of the emitting 
region. 



4 SIGNATURES OF LINE-DRIVEN 
ACCELERATION 

The most striking observation in this flow is that the veloc- 
ity differences between distinct components coincide with 
various doublet velocity splittings (see Fig ^ and Table. 
The most apparent coincidences are seen for the N v dou- 
blet splitting. Both the C iv lines of components c and d 
match very closely. Not only the velocity difference between 
the centroids of the two components is very close to the 
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Table 2. Velocity splitting among distinct components 



Doublet 


splitting 


components 


separation 




(km s^-"-) 




(km s^-"-) 


C IV 


499 


b - c 


495 






e - g 


499 


N V 


962 


c - d 


969 






d - f 


930 






e - h 


924 


O VI 


1650 


a - c 


1660 






c - e 


1625 



N V doublet splitting but also the residual intensities and 
the widths of the components {FWHM ~ 300 km s~^) are 
nearly identical. Very good matchings are also seen for one 
of the lines of e and h and d and f. In addition to this 
we notice that components h (defined by weak albeit broad 
C IVA1550 and O vi absorptions) and j are separated by 
the N V doublet splitting (see Fig. |). For the O vi split- 
ting, coincidences are seen for components a and c, c and 
e (with similar residuals). This is complemented by coin- 
cidences with C IV splitting between b and c but mostly 
between e and g and between component i and diffuse ab- 
sorption component j. 

In a spectral range where crowding of absorption lines is so 
important, chance coincidences are possible. In order to con- 
sider the probability that such coincidences occur by chance, 
we randomly populated the observed velocity range with 10 
components. Using 10® realizations we find a probability of 
6x10-* for the occurrence of, at least two coincidences, each 
with the velocity splittings of the O VI , N v or C iv doublets 
with a matching uncertainty of ±20 km . If we introduce 
the additional condition that at least two coincidences occur 
between three components in O vi splitting (like the one 
seen among a-c and c-e) the probability becomes < 10~^ 
for a wavelength matching uncertainty of ±20 km s"^. 
This low probability for mere velocity coincidence together 
with tight matching of some of the profiles and consis- 
tent covering factors clearly confirms the line-locking sit- 
uation in this system. It must be noted that, although the 
corresponding wavelength range is crowded by intervening 
lines, additional coincidences are apparent between the cor- 
responding O VI absorptions. In addition, the velocity sep- 
aration between components a and the red component of i 
(A„ = 5830 km s"^) is very close to the N v to Lya ve- 
locity splitting. Therefore, all components are connected by 
a complex set of coincidences. Under the assumption that 
line-locking is structuring the flow, it is expected that com- 
ponents with larger covering factors be located closer to the 
source. This would imply in that case that the flow is an 
accelerating wind. 



5 DISCUSSION 

The physical conditions in part of the flow can be inves- 
tigated from the very peculiar component e at Zabs = 2.8606 
(—1350 km in Fig. Absorption lines from H I Ly- 
man series lines as well as excited levels of singly ionized 
species C ll and Si ll are present (see Fig. ^ and from C ll 
to P V. However, none of the Lyman series absorption lines 
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Figure 3. Velocity plot of Lyman series absorption lines pro- 
duced by component e at z^^^^ = 2.8606 toward Q 1511+091. 
The best fitted profile (dotted line) is obtained considering the 
gas covers only 40% of the continuum emission (from the accre- 
tion disk). Note that a weak additional intervening Lyman a line 
is blended with the Lyman /3 line. The effect of partial coverage 
can be appreciated by comparing the profiles of this component 
with that of an intergalactic absorption with similar amount of 
H I seen at 300 km s^^. 



are dark. Lyman-a is not even seen as a distinct compo- 
nent when it should be conspicuous. Similar situation was 
observed by Telfer et al. (1998) toward QSO SBS 1542±531. 
From the flat bottom of the Lyman-/? line, it is apparent that 
the cloud does not cover the background source completely. 
Higher Lyman series lines are well fltted with a covering fac- 
tor of the order of 0.4. It must be noted that the emission 
lines for the corresponding transitions are negligible com- 
pared to the continuum. In addition, such a low covering 
factor can not be accommodated by scattered light (Ogle et 
al. 1999) . Therefore, it can be concluded that the cloud cov- 
ers only about 40% of the central source of continuum whose 
radius is of the order of 10"'^ pc. All neutral hydrogen lines 
can be fitted consistently (see Fig. ^ with /c = 0.4 and 
log 7V(H l) ~ 16.30. A consistent model for all other absorp- 
tion lines is found with this covering factor (see Fig. ^ and 
using the Al ill and Si ill profiles as templates for, respec- 
tively, the singly-ionized and high-ionization species. The 
difference between the two profiles is mainly due to an extra 
component in the red wing that is seen in Si ill and other 
higher ionization lines. Note that although the centroids of 
the lines are the same, the width of the lines increases with 
ionization. 

It is striking to note that with this covering factor, the ratios 
C ll*/C II and Si ll*/Si ll are found to be both of the order 
of two which is the value expected in case of thermal equilib- 
rium. Given the ionization state of the gas (see below), the 
temperature is most probably larger than 10* K, and from 



rn\ nnnn R A « a/twr a s nnn nnn_nnn 



Collimated flow driven by radiative pressure 5 




- ^A111I>,1854 



-200 200 

Relative velocity (kms"') 

Figure 4. Fits of the absorption profiles in component e at 
2abs = 2.8606 toward Q 1511+091. The observed profile of 
Al IIIA1854 is used as a template to fit the singly ionized species. 
The profile are fitted assuming a covering factor of 0.4 and by 
varying the central optical depth preserving the profile function 
given by the template. P V and Si IV profiles are fitted using the 
Si III profile as the reference template. The difference between the 
two profiles is mainly due to an extra component in the red wing 
that is seen in Si in and other higher ionization lines. 



the thermal equilibrium between populations of Si ll levels 
we can derive the electronic density rie > 1.3x 10"^ T^'^ cm~'^ 
(T4 is the temperature in the units of 10^ K). For a Math- 
ews & Ferland (1987) QSO ionizing spectrum, and using 
the constraints from observed ion ratios, the ratio of the 
ionizing photon density to the hydrogen particle density 
(the so-called ionization parameter, U) can be shown to be 
in the range —2.0 < log U < —0.6. The metallicities are 
most probably larger than solar. This is somewhat consis- 
tent with the typical values obtained for the Broad Line 
Region, -1.0 < log f/ < -0.0 and 10^ ^ ^ lO" cm"^, 
although the density is probably smaller in the cloud (see 
below) . Assuming a typical luminosity density at the Lyman 
limit of 10^^ erg s~^ Hz~^, and a spectral index of —1.6 the 
maximum distance of the absorbing gas from the source is 
15.75 Y^l/n?7 kpc. The lower limit on rie derived above sug- 
gests the gas has to be at a distance less than a few kpc 
from the ionizing source. 

From the absence of absorption due to C lll*A1176 origi- 
nating from an excited metastable level, we derive an up- 
per limit on the density of 8x10^ cm~^. This gives a mini- 
mum distance of a few 10~^ pc. Thus the observations are 
consistent with the absorbing gas being inside or close to 
the BLR. For a spherical absorbing cloud the particle den- 
sity is ~ A''(H)/v^i?c, where Rc and Af(H) are the radius 
of the continuum emitting region and the total hydrogen 
column density of the absorbing gas (in the range, 10^^ 



A'^(H)^ 10^" cm~^, from the models) respectively. For the 
lower limit on density we get, Rc ^ 0.04 T^^'^pc. 
The largest blue-shift observed in the C IV profile is much 
smaller than the typical terminal velocities of BAL flows. 
This probably means that the gas could be further accel- 
erated. Indeed, for physical conditions similar to those ob- 
served in different components of the Q 1511+091 outflow, 
the line-radiation force is much larger than gravity even 
in the close proximity of the black-hole (de Kool & Begel- 
man 1995) . Most of the hydrodynamical radiatively driven 
wind models predict the gas to reach terminal velocity at 
a distance of a few pc from the ionizing source (Arav et 
al. 1994). Lack of acceleration in components with velocity 
much lower than the terminal velocity could be accounted 
for if either our line of sight probes a small portion of a 
non-radial flow (Arav 1996) or the extra drag forces are at 
play in reducing the effective acceleration of the gas. In the 
case of Q 1511+091, the absorption lines are seen at redshift 
very close to or slightly higher than the systemic redshift. 
In addition, nearly perfect velocity matching with doublet 
splittings suggests that the direction of the flow is very close 
to our line of sight. 

Disk winds have been suggested as the natural origin for 
ejected material from quasars giving rise to BAL systems 
(de Kool & Begelman 1995; Murray et al. 1995; Proga et al. 
2000). In these models, a cylindrically symmetric disk wind 
is ejected away from the surface of the disk and is radiatively 
accelerated radially once exposed to the strong continuum 
source. Line-driven instabilities could lead to the formation 
of distinct velocity components seen in the flow (Feldmeier et 
al. 1997). Our observations illustrated in Fig. H support the 
basic mechanism at work in these models. It would be very 
interesting to constrain the acceleration over a large time- 
scale for our understanding of the nature of drag forces and 
in general the nature of the flow. 
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Figure 5. Qualitative sketch of the relative position of the absorbing gas compared to the flux emitting region. The fact that some of 
the components of the outflow from Q 1511+091 cover only part of the source of continuum suggests that the absorbing gas is close to the 
continuum emission region. Typically the continuum is believed to originate from a region of size R ~ 150 AU. Reverberation mapping 
studies of AGNs indicate that the broad emission lines are produced by a region of radius at least 100 times that of the continuum 
emitting region. Thus it is most likely that at least part of the observed flow originate close to the inner boundary of the broad emission 
line region. 
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